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We evaluate 0{af,as) corrections in the Higgs boson sector of the CP-conserving 
X ■ MSSM, generahsing the known result in the literature to arbitrary values of tan/3. 

A detailed analysis of the renormalisation in the bottom/scalar bottom sector is per- 
formed. Concerning the lightest MSSM Higgs boson mass, we find relatively small 
corrections for positive fi, while for ^ < the genuine two-loop 0{abCts) corrections 
can amount up to 3 GeV. Different renormalisation schemes are applied and numer- 
ically compared. It is demonstrated that some care has to be taken in choosing an 
appropriate renormalisation prescription in order to avoid artificially large corrections. 
The residual dependence on the renormalisation scale is investigated, and the remain- 
ing theoretical uncertainties from unknown higher-order corrections in this sector are 
discussed for different regions of the MSSM parameter space. 
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1 Introduction 



A crucial prediction of the Minimal Supersymmetric Standard Model (MSSM) [1] is the 
existence of at least one light Higgs boson. The search for this particle is one of the main 
goals at the present and the next generation of colliders. Direct searches at LEP have already 
ruled out a considerable fraction of the MSSM parameter space [2,3], and the forthcoming 
high-energy experiments at the Tevatron, the LHC, and the International Linear Collider 
(ILC) will either discover a light Higgs boson or rule out Supersymmetry (SUSY) as a 
viable theory for physics at the weak scale. Furthermore, if one or more Higgs bosons are 
discovered, bounds on their masses and couplings will be set at the LHC [4-6]. Eventually 
the masses and couplings will be determined with high accuracy at the ILC [7-9]. Thus, a 
precise knowledge of the dependence of masses and mixing angles in the MSSM Higgs sector 
on the relevant supersymmetric parameters is of utmost importance to reliably compare the 
predictions of the MSSM with the (present and future) experimental results. 

The status of the available results for the higher-order contributions to the neutral CP- 
even MSSM Higgs boson masses can be summarised as follows. For the one-loop part, the 
complete result within the MSSM is known [10-13]. The dominant one-loop contribution is 
the 0{at) term due to top and stop loops {at = /if/(47r), ht being the superpotential top 
coupling). Corrections from the bottom/sbottom sector can also give large effects, in partic- 
ular for large values of tan/5, the ratio of the two vacuum expectation values, tan/3 = V2/V1. 
The computation of two-loop corrections is also quite advanced. It has now reached a stage 
such that all the presumably dominant contributions are known. They include the strong 
corrections, usually indicated as 0{atas), and Yukawa corrections, 0{a^), to the dominant 
one-loop 0{at) term, as well as the strong corrections to the bottom/sbottom one-loop 
0{ab) term {ab = hi/ {Air)), i.e. the 0{abas) contribution, derived in the limit tan/3 — ^ 00. 
Presently, the 0{atas) [14-23], 0{af) [14,15,24,25] and the C(a6a,) [26] contributions to 
the self-energies are known for vanishing external momenta. Most recently also the cor- 
rections 0{atab) and 0{al) [27], a "full" two-loop effective potential calculation [28] and 
an evaluation of the leading two-loop momentum dependent effects [29] have become avail- 
able. In the (s)bottom corrections the all-order resummation of the tan /5-enhanced terms, 
0{ab{as tan/5)"), is also performed [30,31]. Reviews with further references can be found in 
Refs. [32-34]. 

The b/b sector has attracted considerable attention in the last years, since its corrections 
to the MSSM Higgs boson sector have been found to be large in certain parts of the MSSM 
parameter space, possibly even exceeding the size of the top/stop corrections. This can 
happen especially for large values of tan f3 and the supersymmetric Higgs mass parameter /x. 
For illustration, we show in Fig. 1 the shift in the lightest CP-even Higgs boson mass, AM^, 
arising from the b/b sector at the one- loop level (all two-loop corrections are omitted here) 
as a function of the bottom-quark mass for large tan/3 and \fi\. The bottom-quark mass 
in this plot is understood to be an effective mass that includes higher-order effects (see the 
discussion in Sect. 3). The figure demonstrates that corrections from the b/b sector can get 
large if the effective bottom mass is bigger than about 3 GeV. 

The possibly large size of the corrections from the b/b sector makes it desirable to in- 
vestigate the corresponding two-loop corrections and thus to analyse the renormalisation in 
this sector. An inconvenient choice could give rise to artificially large corrections, whereas 
a convenient scheme absorbs the dominant contributions into the one-loop result such that 
higher-order corrections remain small. The comparison of different schemes (where no ar- 



1 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 



m^^" [GeV] 

Figure 1: The shift in the hghtest CP-even Higgs-boson mass from the one- loop correc- 
tions in the h/h sector is shown as a function of the (effective) bottom-quark mass for 
H = ±1000 GeV, tan/3 = 50, MgusY = 600 GeV, At = = 500 GeV, Ma = 700 GeV. 



tificially enhanced corrections appear) gives an indication of the possible size of missing 
higher-order terms of 0{abCtl)- 

In this paper we derive the result for the 0{abas) corrections in various renormalisation 
schemes. The relations between the different parameters in these schemes are worked out in 
detail. The absorption of leading higher-order contributions into an effective bottom-quark 
mass is discussed. We perform a numerical analysis of the various schemes and compare our 
results with a previous evaluation of the 0{abas) corrections carried out in the limit where 
tan/3 is infinitely large [26]. We discuss the dependence of our result on the renormalisation 
scale and provide an estimate of the remaining theoretical uncertainties in this sector.^ 

The paper is organised as follows: in Sect. 2 we briefly review the MSSM Higgs boson 
sector, outline the corresponding renormalisation at the two-loop level, and describe the 
evaluation of the diagrams of 0{ab(ys)- Sect. 3 contains a detailed description of the renor- 
malisation of the scalar top and scalar bottom sector, which is explicitly carried out in four 
different renormalisation schemes for the latter. The numerical analysis of the 0{abas) cor- 
rections, the comparison of the different schemes, the investigation of the renormalisation 
scale, and the comparison with the previous result are performed in Sect. 4. The conclusions 
can be found in Sect. 5. 

^This kind of issues have not been addressed in Refs. [28,29]. 
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2 The Higgs sector at higher orders 



We recall that the Higgs sector of the MSSM [35] comprises two neutral CP-even Higgs 
bosons, h and H {rrih < mn), the CP-odd A boson, ^ and two charged Higgs bosons, if^. 
At the tree-level, the masses rrih^tree and m^/ tree can be calculated in terms of Mz, Ma and 
tan (3 from the mass matrix for the neutral CP-even Higgs components (denoted by </>) 

_ / Ml sin^ /5 + M| cos^ /3 -(M^ + M|) sin/3 cos/? \ 

~ V -(Ml + M|) sin /? cos /3 Ml cos2 /3 + M| sin^ /? J 

by diagonalization, 

("^T' ' ) = U,M,Ul U,= ( '^'^ ^^^"V (2) 

V "^Itrec / ^ <P <P <P y -sma cosa J 

with the angle a determined by 



^ M^ + M| „ 
tan 2a = tan 2/5— I -|, -- < a < 0. (3) 



vr 



Ml - M| ' 2 



In the Feynman-diagrammatic (FD) approach, the higher-order corrected Higgs boson 
masses, Mh and Mh, are derived as the poles of the h, if-propagator matrix, i.e. by solving 
the equation 



P 



m 



h,tree 



V 



m 



H,trcc 



-'HH 



^huip^ 



0. 



The renormalised self-energies 




^knip': 

±kh{p'] 



can be expanded according to the one-, two-, . . . loop-order contributions, 

S(/) = S«(/) + S(2)(p2) + -- - . 



(4) 



(5) 



(6) 



The dominant one-loop contributions to the Higgs boson self-energies (and thus to the Higgs 
boson masses) from the b/b sector are of 0{ab) and arise from the Yukawa part of the theory 
(neglecting the gauge couplings) evaluated at p"^ = 0. This has been verified by comparison 
with the full one-loop result from the b/b sector. Hence, the leading two-loop corrections 
from the b/b sector are the 0{as) corrections to those dominant one-loop contributions; 
they are obtained in the same limit, i.e. for zero external momentum and neglecting the 
gauge couplings (the same approximations have been made in Ref. [26]). This approach is 
analogous to the way the leading one- and two-loop contributions in the top/stop sector have 
been obtained, see e.g. Ref. [19]. 

The renormalisation of the Higgs-boson mass matrix for the 0{af,as) corrections under 
consideration follows the description for the 0{atas) terms given in Ref. [19]. Renormalisa- 
tion can be performed by adding the appropriate counterterms. 



,(2) 



(7) 



^Throughout this paper we assume that CV is conserved. 
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where SA4^ denotes the zth-loop counterterm matrix consisting of the counterterms to the 
parameters in the tree-level mass matrix (1). Field renormalisation is not needed for the 
leading 0{af,as) corrections. The renormalised two-loop Higgs boson self-energies with the 
leading contributions of 0{ab(ys) are thus given by 

S(2)(0) = S(2)(0)-W<^5AlfwJ. (8) 



The counterterm matrix in (8) is composed of the counterterms for the A-boson mass 
and for the tadpoles th^n (with s„ = sinOw, Cw = cos6'iy), 

f \^-sm/5cosp cos p J ^ 

+ 7777^ (--os(3{l + sin^ /5) - sin^ ^ \ ^ ^^g) _ ^.^ ^ ^^(,) 
2MzCwSw V -sm^/? cos/5sm^/3y^ " 

+ — n /i 2o^ ■ ^ (sma(5tV + cosa(5t^O- (9) 

2Mzc^Sw V -cos^/^ -(1 + cos2/5)sm/5y ^ " 

The counterterms are determined by the following conditions: 

(i) On-shell renormalisation of the A-boson mass, formulated in the approximation of 
vanishing external momentum, determines the two- loop A- mass counterterm 5Mj^ 
according to 

5Ml^'^ = S!?) (0) . (10) 

(ii) Tadpole renormalisation determines the tadpole counterterms by the requirements 

<' =-<2'. «f =-<f. (11) 

which means that the minimum of the Higgs potential is not shifted. 

The genuine two-loop Feynman diagrams to be evaluated for the Higgs boson self-energies 
and the tadpoles are shown Fig. 2 and Fig. 3. The diagrams with subloop renormalisation 
are depicted in Fig. 4 and Fig. 5. The counterterms for the insertions, where different renor- 
malisation schemes will be investigated, are specified in the next section. 

The diagrams and the corresponding amplitudes have been generated with the package 
FeynArts [36,37]. The further evaluation has been done using the program TwoCalc [38]. 
The resulting expressions are given in terms of the one-loop functions Aq and Bq [39], and 
the two- loop vacuum integrals [40]. 



3 Renormalisation of the quark/squark sector 

Since the two-loop self-energy is considered at 0{a^t,b}(^s) it is sufficient to determine the 
counterterms induced by the strong interaction only. 

The squark-mass terms of the Lagrangian, for a given species of squarks q, can be written 
as the bilinear expression 

C,.^^s = -{qlq^)MJf), (12) 
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Figure 2: Generic two-loop diagrams for the Higgs-boson self-energies (0 = h, H, A; 
hhk.l = 1,2). 
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Figure 4: Generic one-loop diagrams with counterterm insertion for the Higgs-boson self- 
energies (0 = /i, H, A), i, j, k = 1,2). 
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Figure 5: Generic one-loop diagrams with counterterm insertion for the Higgs tadpoles 
(0 = /z,i/, ^,j = l,2). 



with M.q as the squark-mass matrix squared, 

' V - Ml + ml + Mlc,,Q,slJ ' ^ 

where the quantities Ml, M?^, Aq are soft-breaking parameters, and n is the supersymmetric 
Higgs mass parameter. Since we are dealing in this paper with a CP-conserving Higgs sector, 
these parameters are treated as real. As an abbreviation, C2/3 = cos(2/3) is introduced; k is 
defined as k = cot P for up-type squarks and k = tan/3 for down-type squarks. rriq, Qg, and 
are mass, charge, and isospin of the quark q. 

The mass matrix (13) can be diagonalised by a unitary transformation, which in our case 
of real parameters involves a mixing angle 6q, 

=«,- ('>) , u, = C^;- ]';A = ( '°'% . (14) 

q2j KQrJ \Uq21 Uq,J y-smOq cosOql 
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In the (gi, g2)-basis, the squared-mass matrix is diagonal 



V,=UMMl=M' % ] , (15) 



92 



with the eigenvalues m?^ and m?^ given by 



''-^ - 2^^ ' + ^'^^ + + 2^'^ ^2\Ml-Ml + Mlc,,i Tl - 2g;4)[ 

{Ml - Ml + Mlc2f,{T^^ - 2Q,sl)Y + ^^l^^i ' ^^^^ ■ (16) 

The squark-mass matrix can now be expressed in terms of the two mass eigenvalues and the 
mixing angle, yielding 

/ cos^ 9, ml + sin^ 6, m| sin 9, cos % (m| - m|) \ 

<? ~ •/! n r 2 2 \ ■ 2 n 2, 2/i 2 1' ^ ^ 

y sm b'q; COS Og (m|^ - m|J sm % + cos^ Og m^^ y 

3.1 Renormalisation of the top and scalar top sector 

The {t,t) sector contains four independent parameters: the top-quark mass rrit, the stop 
masses m^-^ and m^^, and either the squark mixing angle 6^ or, equivalently, the trilinear 
coupling At. Accordingly, the renormalisation of this sector is performed by introducing 
four counterterms that are determined by four independent renormalisation conditions. 

The following renormalisation conditions are imposed (the procedure is equivalent to that 
of Ref. [41], although there no reference is made to the mixing angle). 

(i) On-shell renormalisation of the top-quark mass yields the top mass counterterm, 

6mt = ^mt [Re E^^ (m^) + Re S^^ (m^) + 2 Re S^^ (m?)] , (18) 

with the scalar coefficients of the unrenormalised top-quark self-energy, St(p), in the 
Lorentz decomposition 

(ii) On-shell renormalisation of the stop masses determines the mass counterterms 

6ml = (^1 ) . = Re ^i,, (m| ) , (20) 

in terms of the diagonal squark self-energies. 

(iii) The counterterm for the mixing angle, Bj;, (entering (17)) is fixed in the following way, 

ReSi^^(mn + ReS,-^^(mn 
Sdt = -^^ttH 2 ' 21) 

2(m? - m? ) ^ ^ 

involving the non-diagonal squark self-energy. (This is a convenient choice for the treat- 
ment of 0{as) corrections. If electroweak contributions were included, a manifestly 
gauge-independent definition would be more appropriate.) 
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In the renormalised vertices with squark and Higgs fields, the counterterm of the trihnear 
couphng At appears. Having aheady specified 69^, the At counterterm cannot be defined 
independently but follows from the relation 



. 2mt{At- fx cot P) 

sm 29i = 2 2 ' 22) 

mi — mi 

tl t2 



yielding 



SA 



1 



1 



mt L2 

1 



sm29t(Sm^^ — Smf^) + cos26'j(m?^ — m?^) 69t 



sin 29f(mr — mr ) Smt 

2mt *^ *i . 



(23) 



This relation is valid at 0{as) since both fi and tan/3 do not receive one-loop contributions 
from the strong interaction. 

3.2 Renormalisation of the bottom and scalar bottom sector 

Because of SU(2)-invariance the soft-breaking parameters for the left-handed up- and down- 
type squarks are identical, and thus the squark masses of a given generation are not inde- 
pendent. The stop and sbottom masses are connected via the relation 

cos^ + sin^ 6'^m?_^ = cos^ 6'^-m| + sin^ 6'im| + ml - m] - cos(2/3) , (24) 

with the entries of the rotation matrix in (14). Since the stop masses have already been 
renormalised on-shell, only one of the sbottom mass counterterms can be determined inde- 
pendently. In the following, the 62 mass is chosen^ as the pole mass yielding the counterterm 
from an on-shell renormalisation condition, i.e. 

Sml = ReE-,Jml) , (25) 

whereas the counterterm for m^^ is determined as a combination of other counterterms, 
according to 

Sml = (cos' 9i6ml + sin' 9i6ml - sin' 9-,6ml - sin 2^,(m| - m|)5^, 

COS V 

-|- sin 26'^(m?^ — m^^59i — 2mt 6mt + 2m{, Sm^j ■ (26) 

Accordingly, the numerical value of m^ does not correspond to the pole mass. The pole 
mass can be obtained from via a finite shift of 0{as) (see e.g. Ref. [42]). 

There are three more parameters with counterterms to be determined: the 6-quark mass 
m^, the mixing angle 9^, and the trilinear coupling Ai,. They are connected via 

. 2mb{Ab- fitanp) 

sm 29i = 2 2 ' 27 

mi — mi 
01 62 



^This choice is possible since (14)-(16) ensure that the 62-field and the foi-ficld do not coincide. 
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which reads in terms of counterterms 



2 cos 26 861 



sin 29' 



5mh 



2mh 5 Ah 



b '2 
rrib mf 

01 



- sin2^r 



5mi 

b2 



62 



2 2 



{21 



b2 



Only two of the three counterterms, (5m;,, 56'^, 6Aiy can be treated as independent, which offers 
a variety of choices. In the following, four different renormalisation schemes, see Tab. 1, will 
be investigated. Two of them are on-shell schemes in the sense that the Higgs self-energies 
do not depend on the renormalisation scale fjP^. 



scheme 


2 

mf 


rrib 


Ab 


Ob 


analogous to t/t sector ("m^ OS") 


on-shell 


on-shell 




on-shell 


DR bottom-quark mass ( "m^ DR" ) 


on-shell 


DR 


DR 




DR mixing angle and Ai, ("Ah, 9^ DR") 


on-shell 




DR 


DR 


on-shell mixing angle and At, ("Afe, 9^ OS") 


on-shell 




on-shell 


on-shell 



Table 1: Summary of the four renormalisation schemes for the bottom quark/squark sector inves- 
tigated below. Blank entries indicate dependent quantities. 

The schemes are described in the following subsections, prior to the discussion of their 
quantitative numerical features in Sect. 4. 

3.2.1 Analogous to the top quark/squark sector 

A straight-forward possibility is to impose renormalisation conditions in analogy to those of 
the top quark/squark sector in Sect. 3.1. 

(i) On-shell renormalisation of the bottom quark mass rrib determines the corresponding 
counterterm as follows, 



6mb = -nib [Re Sb^ {ml) + Re S^^ {ml) + 2 Re (m; 



(29) 



(ii) The counterterm for the sbottom mixing angle 9i is determined in the following way. 



59-. 



ReS^^^(mf^)+ReS^^^(mf; 

2(m? — m? ) 
^01 02 ^ 



(30) 



The dependent counterterm (5m? for the hi mass is then fully specified by (26). Moreover, 
Ab is treated here as a dependent quantity; the corresponding counterterm 5Ab follows from 
the relation (28), yielding in combination with (26) the expression 



5Ab 



nib 



- tan 9i5ml^ + (m?^ - m|)5^^ - (5mt 



2mj, 



(m? — m?_^ ) sin 26*^ — 2 tan 6'^m;, 



-|- tan 6'^!^ cos 9i5m^^ + sin 9^5m^^ — sm29^{m^^ — m^^)59i — 2mt6mt 
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(31) 



While formally the renormalisation described in this section is the same as in the top/stop 
sector, there are nevertheless important differences. The top-quark pole mass can be directly 
extracted from experiment and, due to its large numerical value as compared to other quark 
masses and the fact that the present experimental error is much larger than the QCD scale, 
it can be used as input for theory predictions in a well-defined way. For the mass of the 
bottom quark, on the other hand, problems related to non-perturbative effects are much more 
severe. Therefore the parameter extracted from the comparison of theory and experiment [43] 
is not the bottom pole mass. Usually the value of the bottom mass is given in the MS 
renormalisation scheme, with the renormalisation scale /i^^ chosen as the bottom-quark 



Another important difference to the top/stop sector is the replacement of cot (3 — ^ tan (3. 
As will be discussed in more detail below, very large effects can occur in this scheme for 
large values of /i and tan /5. 



3.2.2 DR bottom-quark mass 

Potential problems with the bottom pole mass can be avoided by adopting a renormalisation 
scheme with a running bottom-quark mass. In the context of the MSSM it seems appropriate 
to use the DR scheme [44] and to include the SUSY contributions at 0{as) into the running. 
We therefore choose a scheme where m;, and Ai, are both renormalised in the DR scheme. 
The following renormalisation conditions are imposed for the independent quantities. 

(i) The 6-quark mass is defined in the DR scheme, which determines the mass counterterm 
by the expression 



where S'^'^ means replacing the one- and two-point integrals A and Bq in the quark 
self-energies by their divergent parts in the following way. 



with A = 2/e - 7 log47r, and D = A - e. 

(ii) Besides m;,, also the trilinear coupling A;, is defined within the DR scheme. Using (31) 
and inserting the self-energies yields the counterterm 



mass, i.e. mf^{mf^) [43]. 



5m, = \m, [Re T.f^{ml) + Re T.f;{ml) + 2 Re S^f (m^)] , 



(32) 



^("^)|div = "^^A , 
5o(J9^ml,m2)|div = A , 



(33) 





- ml (Re S^^(m,2) + Re T.'^^{ml) + 2 Re T.f^{ml)) 




(Re T.f^{ml) + Re S^-(m^) + 2 Re T.f^{ml)) . 



(34) 
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The counterterms for the mixing angle, 59i, and the hi mass, 5m? , are dependent quantities 
and can be determined as combinations of the independent counterterms, invoking (26) 
and (28), 

rrLhSAh + tan6'^(5m|_^ + 6mi,(^- — (m?^ — mj^^) sin 26'^ — 2tan^^mb^ 



S9z 



b 2 2 

m~ — m~ 

01 02 



-2m;, 

— tan 9i ^cos^ 9i5m^^ + sin^ 9i5m^^ — sin 26'j(m?^ — m^^_^)59i — 2mtSmt^ , (35) 

(5m?^ = tan^ 9i5rn^^ + 2 tan^^mb(5/lfe + 2^ — sin^ 9i{m^^^^ — m?_^) + (1 — tan^ 9i)mbj6mij 

+ (1 - tan^ 61^) (cos^ 9i5m\^ + sin^ 9i5m\^ - sin 2%(m| - m]^)59i - 2mt5mi^ . (36) 



The renormahsed quantities in this scheme depend on the DR renormahsation scale fi^^. 
If not stated differently, in all numerical results given in this paper the DR scale refers to 
the top-quark mass, i.e. = m^. 

In order to determine the value of m|;^^'^^^^(/i°^) from the value m^^(/i^^) that is 
extracted from the experimental data one has to note that by definition jjj^^'^^^^ contains 
all MSSM contributions at 0{as), while m^^ contains only the 0{as) SM correction, i.e. the 
gluon-exchange contribution. Furthermore, a finite shift arises from the transition between 
the MS and the DR scheme. As input value for m^^(M^) we use in this paper m^^{Mz) = 
2.94 GeV [45]. _ _ 

The expression for m^^'^^^^(/i°^) is most easily derived by formally relating 77^^^'^^^^ 
to the bottom pole mass first and then expressing the bottom pole mass in terms of the MS 
mass (the large non-perturbative contributions affecting the bottom pole mass drop out in 
the relation of mf^'^^^^ to mp). Using the equality m^^ + Smf^ = mp'^^'^^ + ^mf^'^^^^ 
and the expressions for the on-shell counterterm and the DR counterterm given in (29) and 
(32), respectively, one finds 

^DR.MSSM^^i5K) ^ ^os + l^,(s«^^(W) + S«^^(m,^)) + m, s£(m^) . (37) 

Here the T,^"^ are the UV-finite parts of the self-energy coefficients in (29). They depend 
on the DR scale /i^^ and are evaluated for on-shell momenta, = m^. Inserting m^^ = 



m^^{Mz)b'^'^\ where 



Lshift 



MS\2 



— - - In 



TT V3 m: 



2 



(3^ 



one finds the desired expression for m^^ 



DR,MSSM(^DR^ = mp(Mz)6^^'f* + ^m, (s«^-(m,2) + ^I^M) + m, S£(mg) . (39) 



3.2.3 DR mixing angle and Ab 

A further possibility is to impose renormahsation conditions for the mixing angle 9^ and for 
Ah, and to treat the counterterm of the 6-quark mass as a dependent quantity determined 
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as a combination of the other counterterms using the relation (28). The renormalisation 
conditions in this case read explicitly: 



(i) 6Ab is determined in the DR scheme as in the previous case by the expression (34). 

(ii) The mixing angle 6^, defined in the DR scheme, is renormalised by the counterterm 



ReSf^(m? ) +ReSf^(m? , 

012 01 012 02' 



69-, = "-^^^^ ^ — — . (40) 



2(mf — mf , 

^ bi b2' 



The counterterm for the fe-quark mass, ^m^, can be obtained using (28) and the con- 
straint (26). It is given by the following quantity (which is well-behaved for — *■ 0), 



6mh 



tan 9i ^— (5m^^ + cos^ 9f6mj^ + sin^ 9i6rn?^^ — sin 29i{m^^^ — m^^)69i — 2mt 6mt^ 



2 2 

mi — m 



- mbMfc + (m?^ - m|)(5^^J ^— ^ sin 2^^ - 2mbtan^^J . (41) 

The numerical value of m^ in this scheme is obtained from (39) and the (finite) difference of 
the counterterms given in (41) and (32). 

Finally, (26) yields also the counterterm for the dependent squark mass, 5m? , with the 
specification (41) for the 6- mass counterterm. 

3.2.4 On-shell mixing angle and Ai, 

In Ref. [26] a renormalisation condition was imposed on the Abib2 vertex in order to avoid 
an explicit dependence on the renormalisation scale For the purpose of comparing our 
results with those of Ref. [26] we include such a renormalisation scheme in our discussion. 
While in Ref. [26] the limit tan/3 oo has been used to derive all the renormalisation 
conditions and counterterms, we have derived the relevant quantities for arbitrary values of 
tan/3. We call this scheme "on-shell" (as in Ref. [26]), although the vertex is taken at an 
off-shell value of the A-boson momentum. 

Similarly to the previous scheme, the counterterm for the 6-quark mass is derived as a 
linear combination of other counterterms by means of (28). The independent renormalisation 
conditions can be formulated as follows. 

(i) The counterterm for the mixing angle 9^ is defined by 

Re Er (m? ) + Re Sr (m? ) 
69, = ^^^^ , (42) 

^ bi 02 

as in the scheme "analogous to the top quark/squark sector". 

(ii) Ab is determinded by imposing the condition 

A(0,m?^,m[) + A(0,m?^,m?J = , (43) 

with A(p^,p? ,p?^) as the renormalised three-point Ahih2 vertex function. 
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A 



■o 



le 



mfetan (3 6 A}, 



(44) 



In the large-tan/? limit, this requirement reproduces the condition apphed in Ref. [26]. 

Condition (ii) can be formulated as an equation determining the counterterm for A^ 
in the following way, 



5A 



b = 



+ 



. Mw sin 9 



e mi, tan 
fi + A{, tan P 
nib tan j3 



f (A(0,m^^,m^^)+A(0,m4,m4 
. Mw sin 9w 



fi + AhtanP 
2 tan/? 



62' 



mb(/i + tan /3) \ , / 2 2 n./j 

+ l^^p + + K, - "^te)^^'^ 

— tan6'^^5m?_^ — cos^ O^Sm'j^ — sin^ O^Sm"^^ + sin26'f(m?^ — m^SO^ 
+ 2mt5mt 
where the Z factors are defined as 



/i(tan/3 H ;-) + 2mbtan6'r 

^ tanp^ 



(45) 



bibi 



Sr (m? ) - Sr (m? ) 

2 2 



(46) 



61 



b2 



Again, the dependent counterterm for the 6-quark mass is determined by (28) and the 
constraint (26), but now inserting the above specification (45) for 6Ab, yielding 



Smh 



. Mw sin 9w 

e tan [3 



A(0,m?^,mf^) +A(0, m?^,mp) + 



mf,(/i + Afotan/?) 
2tan/3 



fete ^ 



+ (m?^ - mldOi + tan%(^- sm2ei{m\ - mldOi - Smf^ + cos^ 9iSm~^ 



+ sin^ 6i5m^^ — 2mt6mt 



1 T -1 

/i(tan/3 H -) + 2mbtan6'r 

. ^ tanp^ 



(47) 
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The numerical value of in this scheme is obtained from (39) and the (finite) difference of 
the counterterms given in (47) and (32). 

With the specification of Sruh in (47), also the 6i-mass counterterm 6m? in the general 
relation (26) is fully determined. 

3.3 Resummation in the b/b sector 

The relation between the bottom-quark mass and the Yukawa coupling hf,, which in lowest 
order reads m^ = hbVi/^/2, receives radiative corrections proportional to hbV2 = /ife tan/3 t^i. 
Thus, large tan /3-enhanced contributions can occur, which need to be properly taken into 
account. As shown in Refs. [30,31] the leading terms of 0{ai,{astan (3)"') can be resummed 
by using an appropriate effective bottom Yukawa coupling. 

Accordingly, an effective bottom-quark mass is obtained by extracting the UV-finite 
tan /^-enhanced term Ami, from (39) (which enters through E^^) and writing it as 1/(1 -|- 
Amj,) into the denominator. In this way the leading powers of (ctstan/?)" are correctly 
resummed [30,31]. This yields 

mMSSM I5R ^r(Mz)fc^^^^* + >6(sfi°(ma + S£(W)) +m,S£(m2) 

' ^ = TTa^, ' ^^^^ 

where S^g = + Am;, denotes the non-enhanced remainder of the scalar 6-quark self- 
energy at 0{as), and W^^^^ is given in (38). The tan /3-enhanced scalar part of the 6-quark 
self-energy, Am^, is given at 0{as) by^ 

2 

Amb = —as tan P fxmg I (ml , ml , ml), (49) 



with 



m? m? log(m? /m? ) + m? m?log(mj /m?) + m^ml log(m?/m? ' 

2 2 2\ _ 6i b2 ^ ^2' bi'^ hx 9 '^^ hi' 9' ' 9 b2 '^^ 9' h2' 



limr , mr , m- 



\ a'\ a h^'\ bi' 



(50) 



and Amb > for yU > 0. 

In the "mfo DR" scheme we use the effective bottom-quark mass as given in (48) ev- 
erywhere instead of the DR bottom quark mass (in particular, we use this bottom mass 
in the sbottom-mass matrix squared, (13), from which the sbottom mass eigenvalues are 
determined). The numerical values of the bottom-quark mass in the other renormalisation 
schemes can be obtained from (48) as explained above, and from (51) below. 

We incorporate the effective bottom-quark mass of (48) (or the correspondingly shifted 
value in the other renormalisation schemes) into our one-loop results for the renormalised 
Higgs boson self-energies, which determine the Higgs boson masses at one-loop order accord- 
ing to (4)-(6). In this way the leading effects of 0{abOis) are absorbed into the one-loop 
result. We refer to the genuine two- loop contributions, which go beyond this improved 
one-loop result, as "subleading 0{abOis) corrections" in the following. 

^There are also corrections of 0{at) to Amb that can be resummed [31]. These effects usuahy amount up 
to 5-10% of the 0[as) corrections. Since in this paper we are interested only in the 0[abas) contributions 
to the MSSM Higgs sector, these corrections have been neglected. Further corrections from subleading 
resummation terms can be found in Ref. [46]. 
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4 Numerical results 



4.1 Evaluation 

If not mentioned explicitly in the text the default set of parameters shown in Tab. 2 is used. 
Large values of tan/3 and are chosen in order to illustrate possibly large effects in the h/h 
sector. 



SM parameters: 






mt = 174.3 GeV, mf^{Mz) 


= 2.94 GeV, 


Mz = 91.1875 GeV, Mw = 


80.426 GeV, Gp = 1.16639 10"^ 


parameters of the Higgs sector: 




Ma = 120 GeV 




tan/3 = 50 = -1000 GeV 


soft-breaking parameters: 


for the gauginos: 




for the sfermions: 


5 sin^ 9w , ^ 
Mi= ^ ^ Ms 
3 COS"' 9w 




Ml = Ml - = 1000 GeV with z = 1, 2, 3 


M2 = 100 GeV 




Mj^ = 1000 GeV with f = u, c, t, d, s, b, e, fi, r 


M3 = 1000 GeV 




A{u,c,t} = A{d,s,b} = ^{e,/.,r} = 2000 GeV 



Table 2: Set of default input parameters. 



We will mostly discuss the case of negative /i, since according to Eqs. (48)-(50) this sign 
of /i leads to a negative Arrib and therefore to an increase of the effective bottom-quark mass. 
This gives rise to an enhancement of the corrections from the b/b sector, see Fig. 1. While 
the negative sign of /i is disfavoured from the comparison of the MSSM prediction [47,48] 
with the experimental data on the anomalous magnetic moment of the muon [49], it would 
seem premature at this stage to completely disregard this possibility. For /i > 0, on the other 
hand, the corrections to the Higgs-boson masses from the b/b sector will normally not exceed 
the GeV level if the result is expressed in terms of an appropriately chosen running bottom- 
quark mass (see Fig. 1). It should be noted, however, that the prospective experimental 
accuracy on Mh at the LHC and the ILC will be significantly below the GeV level, so that 
the inclusion of non-enhanced two-loop corrections will be necessary in order to achieve the 
same level of precision for the theoretical prediction (see the discussion below). 

For the calculation of the Higgs boson masses presented below the complete one-loop 
self-energies have been used, with tan j3 renormalised in the DR scheme [50-52] and with 
the Z boson mass on-shell. At the two- loop level, besides the 0{abas) corrections also 
the contributions 0{atas) using the one-loop sub-renormalisation of Sect. 3.1 have been 
included. For simplicity we have neglected the 0{a^) terms [25]. For the 0{atas) corrections 
the top pole mass, = 174.3 GeV, has been used. The inclusion of all known corrections 
and the new experimental top quark mass value of mt = 178.0 GeV [53] in our analysis 
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would yield an increase in of (9(8 GeV) [32]. Therefore the mass values given in our 
numerical analysis should not be viewed as predictions of M/j; they are rather illustrations of 
the ttcj-corrections to the bottom Yukawa contributions at the two-loop level. (It should be 
noted that the chosen parameters are such that they are not in conflict with the experimental 
lower bounds on Mh [2,3].) 

4.2 Comparison of the different renormalisation schemes 

In order to compare the different renormalisation schemes, the parameters entering the one- 
loop result have to be transformed according to the different renormalisation prescriptions. 
As our default for which the input parameters are fixed we have chosen the "m^ DR" scheme, 
where mi, and A}, are defined as DR parameters. As explained in Sect. 3, the parameters are 
converted to a different renormalisation scheme RS (with counterterms with the help 

of the following transformations, 

= mp-(5mf I finite, (51) 

= ^P-5^f|finite, (52) 

and analogously for the other parameters. If not stated differently, the DR scale has always 
been chosen as /i^^ = rrit. As an example, in Tab. 3 and Tab. 4 numerical values for the 
bottom quark mass, Ai, and the sbottom masses in the different schemes (see Tab. 1), are 
shown for tan /3 = 30 and tan /5 = 50 and using the default values given in Tab. 2 otherwise. 



scheme 


nib [GeV] 


Ab [GeV] 


^b. [GeV] 


mfo DR 


3.79 


2000.00 


1059.95 


Ab, Oi DR 


3.04 


2000.00 


1039.50 


Ab, Oi OS 


2.99 


2332.81 


1039.04 


mfo OS 


3.77 


-4284.56 


1039.25 



Table 3: Values of the bottom quark mass, Ab and in the different schemes for tan /3 = 30 
and /i = —1000 GeV. The value of m^^, which is renormalised on-shell (see (25)), is the same 
in all four schemes, m^^ = 938.44 GeV. 

The values given in Tab. 3 and Tab. 4 indicate that the "m^ OS" scheme leads to huge 
corrections in Ab that invalidate the applicability of this scheme. The other schemes give 
rise to only moderate shifts in the parameters. 

The reason for the problematic behaviour of the "m^ OS" scheme is easy to understand. 
The renormalisation condition in the "m;, OS" scheme is a condition on the sbottom mixing 
angle 9i and thus on the combination [Ab — /itan/5) (see (28)). In parameter regions where 
/itan/3 is much larger than Ab-, the counterterm 5Ab receives a very large finite shift when 
calculated from the counterterm 59i. More specifically, 5Ab as given in (31) contains the 
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scheme 


rub [GeV] 


Ah fGeVl 


mr fGeVl 


m, TTR 


o.oz 


9000 00 




At, Oi DR 


5.26 


2000.00 


1117.93 


A, OS 


5.24 


2219.40 


1118.02 


mfo OS 


4.93 


6508.12 


1122.04 



Table 4: Values of the bottom quark mass, A^ and m^^ in the different schemes for tan = 50 
and /i = —1000 GeV. The value of m^^ is the same in all four schemes, m^^ = 836.48 GeV. 

contribution 

6Ab = — 
nib 

= — [-6mb{Ab - fi tan + ...] , (53) 
nib 

that can give rise to very large corrections to Ab. This problem is avoided in the other 
renormalisation schemes introduced in Tab. 1, where the renormalisation condition is applied 
directly to Ab, rather than deriving SAb from the renormalisation of the mixing angle. 

We now turn to the numerical comparison of the different renormalisation schemes. As 
discussed above, the tan /3-enhanced terms of 0{abas) entering via Anib have been absorbed 
into the one- loop result. The meaning of the various curves in the following figures is specified 
as (see also Tab. 1): 

• dashes with dots (black): 0{atas) with m^^'^^^^, results without subleading two-loop 
0{abas) terms 

• dot-dash (light blue): "mfc OS" scheme for subleading two-loop 0{ab<ys) terms 

• solid (red): "mf, DR" scheme for subleading two-loop 0{abas) terms 

• dotted (dark blue): ''Ab, 9^ DR" scheme for subleading two-loop 0{abas) terms 

• dashes with stars (green): ''Ab, 9^ OS" scheme for subleading two-loop 0{abOis) terms 

We start our analysis of the different renormalisation schemes by comparing the results 
for Mh and Mh as a function of tan/3 in Fig. 6. The other parameters are as given in 
Tab. 2. As expected from the discussion of Tab. 3 and Tab. 4, the "nib OS" scheme gives 
rise to artificially large corrections and shows very large deviations from the other schemes 
for intermediate and large values of tan (5. This behaviour is even more pronounced for Mh 
than for M^, as can be seen in the lower plot of Fig. 6. These extremely large corrections are 
a consequence of the large contributions to the counterterm of the parameter Ab (see (53)). 
The Higgs self-energy contribution from virtual sbottoms contains a term proportional to 
A^. Using as input a value for Ab according to (52), very large contributions proportional 
to {SAb)'^ are introduced. These corrections are more pronounced in T^hh, where they enter 
like (cos a AbY, than in T^hh, where they enter like (sin a AbY (1^1 <^ 1 in our analysis). The 



5mb 

2 TUb 



^ bi 



"^62 ) sin 29^ 
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Figure 6: tan /^-dependence of Mh and Mh for Ma = 120 GeV and fi negative. 



unacceptably large contributions to SAf, in the "mfe OS" scheme invahdate a perturbative 
treatment in this scheme. We therefore discard this scheme in the following and focus our 
discussion on the other three schemes defined in Tab. 1. 

The other schemes all give similar and numerically well-behaved results, where Mh starts 
to decrease rapidly with tan /? for tan > 40. Negative mass squares are reached at tan (3 ~ 
53. The main effect comes from the leading contributions of 0{af,as) that enter via the 
resummation of Arrib, see (48). The decrease with increasing tan/5 is mainly due to the 
dependence of Am^ ~ /xtan/5 in (49). The subleading 0{abas) corrections, which arise 
from the genuine two-loop diagrams, are of 0{1 GeV). The differences between the three 
renormalisation schemes are of similar size. For this particular parameter choice the "Af,, 
6i DR" scheme enhances Mh, whereas the other two schemes decrease Mh compared to the 
case where the genuine two-loop corrections are omitted. 
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Figure 7: tan /^-dependence of Mh for Ma = 700 GeV and /i negative. 



In Fig. 7 we show Mh as a function of tan/5 for the same parameters as in Fig. 6, but 
with Ma = 700 GeV. This results in general in larger Mh values, but the general behaviour 
as a function of tan/3 is the same as for Ma = 120 GeV; Mh drops steeply for large tan/3 
values. In all three schemes the subleading terms increase Mh by a few GeV, depending on 
tan j3. 

As discussed above, large corrections from the b/b sector are only expected for negative 
values of /i. In Fig. 8 we show the results for Mh as a function of tan/3 with positive fi 
and Ma = 120 and 700 GeV, respectively. The other parameters are given in Tab. 2. The 
positive sign of /i results in a positive Anib and thus a smaller numerical value of m^^''^^^^. 
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Figure 8: tan /3-dependence of Mh for Ma = 120 GeV (upper pfot) and Ma = 700 GeV 
(fower pfot) for positive /i. 
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As expected^, the variation of with tan/3 is much smaller than for negative /i. Both, 
the leading corrections, i.e. the tan/3 enhanced terms of 0{ai,as), as well as the subleading 
corrections are at the level of (9(100 MeV). The "m^ DR" scheme does not show any visible 
corrections beyond the resummed contributions. This leads to the conclusion that for positive 
fi the corrections beyond the one- loop level coming from the h/h sector are sufficiently well 
under control. However, in view of the fact that the anticipated ILC accuracy on [7-9] 
and the parametric uncertainty of the theory prediction from the ILC measurement of the 
top-quark mass [54,55] will both be about 100 MeV, ultimately the aim will be to reduce the 
theoretical uncertainties from unknown higher-order corrections to at least this level. This 
will require the inclusion of all two-loop corrections (and a significant part of corrections 
beyond two-loop order). For the further analysis in this paper we focus on negative values 
of /i. 

The variation of with (for < 0) for tan/3 = 50 is shown in Fig. 9. As can be 
expected from (49) the corrections at O^ahCXs) increase with increasing Typically the 
genuine two-loop contributions are of 0(1 GeV). For large Ma all the schemes lead to an 
increase of Mh-, whereas for small Ma both negative and positive shifts can occur. Differences 
in the M^ predictions induced by the different renormalisation schemes are below the GeV 
level for large Ma- 

In Fig. 10 the dependence of M^ on Ma is shown for the different renormalisation schemes, 
with the other default parameters from Tab. 2. For Ma ^ 200 GeV the subleading terms 
of all three schemes enhance Mh by 0{1 GeV). A decrease only occurs for small values of 
Ma, depending on the scheme. The differences in the Mh prediction resulting from the use 
of different renormalisation schemes decrease for Ma ^ 200 GeV to (9(0.1 GeV). 

In Fig. 11 it can be seen that the behaviour of the corrections strongly depends on the 
choice of rrig. The figure shows Mh as a function of rUg for // = —1000 GeV, tan/3 = 50 
and Ma = 700 GeV. For nig ^ 1000 GeV all schemes lead to an increase of Mh from 
the subleading 0{a}yas) corrections. For rrig ^ 1500 GeV, on the other hand, all schemes 
lead to a decrease, where the size of the individual corrections also strongly varies with rrig. 
Accordingly, the relative size of the corrections in the different schemes also varies with rrig. 
Corrections up to about 3 GeV are possible. The differences between the three schemes are 
of 0{2 GeV) for large rrig. It should be noted that the effects of the higher-order corrections 
to Mh do not decouple with large rrig. The corrections at 0{atas) [19] as well as 0{ahOis) 
grow logarithmically in the renormalisation schemes that we have adopted. 

The above analysis of the three schemes "mt DR", "A;,, O-^ DR", and "At, 6i OS" in vari- 
ous parameter regions yields numerically well-behaved and physically meaningful results. As 
there is no clear preference for one of the schemes on physical grounds, the difference between 
the results obtained in the three schemes can be interpreted as an indication of the possible 
size of missing higher-order corrections. The size of the individual corrections and also the 
differences between the renormalisation schemes sensitively depend on the input parameters. 
Typically we find that the genuine two-loop corrections in the h/h sector yield a shift in Mh 
of 0(1 GeV). The differences between the three schemes are usually somewhat smaller. 

^See also the discussion in Ref. [26], where the opposite sign convention for is used. 
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Figure 9: /x-dependence of for Ma = 120 GeV (upper plot) and Ma 
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Figure 10: M^i-dependence of Mh for tan/3 = 50 and /i negative. 
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Figure 11: m^-dependence of Mh for Ma = 700 GeV, tan/3 = 50 and /i negative. 
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4.3 Numerical analysis of the renormalisation scale dependence 



While in the previous section we compared the results of different renormalisation scheme, we 
now focus on the "mf, DR" scheme and investigate the effect of varying the renormalisation 
scale of the 0{abas) result obtained in this scheme. We vary the scale within the interval 
< ijP^ < 2mf, resulting in a shift which is formally of 0{abal). The results are shown 
as a function of m-g for tan/3 = 50 in Fig. 12 for Ma = 120 GeV and Ma = 700 GeV. 
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Figure 12: /jP^ dependence of M^ as a function of nig for Ma = 120 GeV (upper plot) and 
Ma = 700 GeV (lower plot) for /i = —1000 GeV, tan/3 = 50. The black area corresponds to 
the 0{atas) result including resummation, i.e. the result without the subleading two-loop 
0{abas) terms. 
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The ijP^ variation of the leading contribution (the 0{atas) result including resummation) 
is shown as the dark shaded (black) band. The results including the subleading corrections 
in the "m;, DR" scheme are shown as a light shaded (red) band. It can be seen that the 
variation with is strongly reduced by the inclusion of the subleading contributions. The 
variation with f/^^ within the "mfe DR" scheme is tiny for rrig ^ 500 GeV, and reaches 
±2 GeV for large rrig values. Thus, the /i^^ variation causes a similar shift in Mh as the 
comparison between the three renormalisation schemes discussed above. 

We have also analysed the variation with fi^^ in the case yU > 0, which is not shown here. 
As for negative /i, the variation with f/^^ is of the same order as the differences between the 
three renormalisation schemes, see Fig. 8. Therefore, for > the unknown higher-order 
corrections to M^, from the b/b sector can be estimated to be of 0(100 MeV). 

4.4 Comparison with existing calculations 

Finally we compare our result with the existing calculation of the 0{abas) corrections pre- 
sented in Ref. [26]. The renormalisation employed there consists of an on-shell renormalisa- 
tion of the two scalar bottom masses and the on-shell condition for Ab shown in Sect. 3.2.4. 
We denote it as "m^, Ab OS" renormalisation. Thus, the differences between our ''Ab, 9^ OS" 
and the "mj^, Ab OS" renormalisation are the different treatment of the m^^ renormalisation, 
as well as the treatment of tan j3. We kept tan as a free parameter, whereas in Ref. [26] it 
was set to infinity in the subleading 0{ab<ys) corrections. In Ref. [26] the shift of the sbottom 
masses due to the SU(2)-invariance was taken into account in the numerical evaluation of 
the sbottom masses following the prescription in Ref. [56] (see also Ref. [42]). 

Our result for in the "A^, $1 OS" scheme is compared with the result of Ref. [26] 
in Fig. 13. For the implementation of the latter ("m^, Ab OS" scheme) the Fortran code 
of Ref. [26] for the numerical evaluation of the 0{asab) corrections to the Higgs-boson self- 
energies has been used [57]. Thereby the input values were determined according to (51) 
and (52). Using these input values for Ab and rrib the sbottom masses were calculated 
taking the sbottom mass shift into account [56]. Mh. is shown as function of rrig for < 0, 
tan/3 = 50, and = 700 GeV. Our result in the "A^, 6^ OS" scheme is shown as the 
dash-star (green) curve, while the result of Ref. [26] ("?ti^, Ab OS" scheme) is given by the 
fine-dotted (pink) curve. The leading contribution in the two schemes, i.e. the 0{atas) 
result including resummation, is also shown: the light-dot-dashed (orange) curve shows the 
0{atas) result using the ''Ab, 9i OS" renormalised parameters; the corresponding result for 
the "mi, Ab OS" renormalised parameters is shown as the light-dotted (gray) curve. 

Fig. 13 shows that the 0{atas) results in the two schemes differ from each other by up 
to 2 GeV for large rrig. The inclusion of the subleading two- loop corrections reduces this 
difference significantly. Our result in the "Ab, 9^ OS" scheme agrees with the result of 
Ref. [26] to better than 0.5 GeV. 

In Fig. 14 we compare our result in each of the three schemes discussed above, i.e. the 
"Ab, Oi OS", the "nib DR" and the "Ab, 6^ DR" schemes, with the result of Ref. [26]. The 
difference AMh between our result and the result of Ref. [26] is shown for each of the three 
schemes as a function of tan/3 for rrig = 1500 GeV, /i = —1000 GeV, and Ma = 700 GeV. 
The differences stay below 1 GeV for tan (3 ^ 50, where our result in the "Ab, DR" scheme 
shows the biggest deviation from the result of Ref. [26] , while as expected, the difference is 
smallest for the "Ab, 9i OS" scheme. For tan /? > 50 large deviations occur because of the 
sharp decrease of M^ in this region (see e.g. Fig. 7). 
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Figure 13: Comparison of our 0{ah<ys) result for Mh in the "A;,, 6^ OS" scheme and the 



result of Ref. [26] 



■ m 



Ai, OS" scheme) as a function of rrig. The 0{atas) results in the 



two schemes, where the subleading 0{ai)as) corrections are omitted (using the appropriate 
renormalised parameters), are also shown. 



<! 



20 



^ = -1000 GeV 
m. = 1500 GeV 
Ma = 700 GeV 




0{abas) ■ m^DR 
Oiatas) : A^, % OS 
: Ah, Oi DR 



25 



30 



35 40 

tan 13 



45 



50 55 



Figure 14: Comparison of our 0(ab«s) result in three different renormalisation schemes and 
the result of Ref. [26]. The three curves for AMh show the difference between our result in 
each of the three schemes and the result of Ref. [26] as a function of tan j3. 
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5 Conclusions 



We have obtained results for the two-loop 0{ab(ys) corrections to the neutral CP-even Higgs- 
boson masses in the MSSM within different renormalisation schemes. The leading tan/3- 
enhanced contributions of the b/b sector can be incorporated into an appropriately chosen 
bottom Yukawa coupling, for which we use the bottom-quark mass in the DR scheme with 
a resummation of the leading contributions. We have analysed in detail the impact of the 
genuine (subleading) 0{ab(ys) two-loop corrections in different parts of the MSSM parameter 
space and we have compared the results obtained in the different schemes. 

We have shown that an on-shell scheme that is frequently used in the t/t sector leads 
to numerically unstable results if it is applied in the b/b sector. The origin of the huge 
corrections in this scheme was traced to the fact that it involves a renormalisation condition 
for the sbottom mixing angle, O^, rather than for the trilinear coupling, Af,. 

The other three schemes that we have analysed yield numerically well-behaved and phys- 
ically meaningful results. For /i > the effect of the genuine 0{af,as) two-loop corrections is 
rather small, typically of 0(100 MeV). Corrections at this level will nevertheless be relevant 
in view of the prospective accuracy of measurements in the Higgs sector and of the top-quark 
mass at the ILC. For < the effective bottom Yukawa coupling increases, leading to an 
enhancement of the effects from the b/b sector. While the constraints arising from the mea- 
surement of the anomalous magnetic moment of the muon favour the positive sign of fi, it 
seems premature at the present stage to discard the parameter region with /i < 0. For large 
values of tan (3 and rrig and large negative values of /i we find that the genuine 0{abas) cor- 
rections can amount up to 3 GeV. We have compared our result for the 0{abas) corrections 
with the existing result in the literature, which was obtained in the limit of tan /3 — >■ oo, and 
found good agreement. 

The comparison of the results in the different schemes that we have analysed and the 
investigation of the renormalisation scale dependence give an indication of the possible size 
of missing higher-order corrections in the b/b sector. For /i > the higher-order corrections 
from the b/b sector (beyond 0{abas)) appear to be sufficiently well under control even in 
view of the prospective ILC accuracy. This applies especially to the "m;, DR" scheme, where 
the corrections beyond the improved one-loop result have been found to be particularly 
small. For /i < 0, on the other hand, sizable higher-order corrections from the b/b sector are 
possible. The size of the individual corrections and also the difference between the analysed 
schemes varies significantly with the relevant parameters, /i, tan/3, rrig and Ma- We estimate 
the uncertainty from missing higher-order corrections in the b/b sector to be about 2 GeV 
in this region of parameter space. 

The results obtained will be implemented into the Fortran code FeynHiggs [58,59]. The 
evaluation of the results within the three schemes will allow to obtain an estimate of the size 
of the missing higher-order corrections as a function of the chosen input parameters. 
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Appendix: Counterterms of the quark/squark sector 

In section 3 the counterterms have been given using the definitions (14) and (16) for the 
sfermion masses and mixing angles. In this appendix the counterterms are given in a more 
general way allowing to use also other definitions for the sfermion masses and mixing angles. 
Introducing a counterterm for the mixing angle needs a certain choice of definitions of the 
sfermion masses and mixing angles. Instead of using an explicit mixing angle counterterm 



the counterterm SYg is introduced as 



= {U,6MM)i2 = {Ug5MM)2i 



(54) 



where the counterterm mass matrix 6A4q contains the counterterms of the parameters ap- 
pearing in (13). With the definitions (14) and (16) 6Yq is related to the mixing angle 
counterterm as follows 



SY, 



(m| -ml)69q 



(55) 



• Top quark/squark sector: 

The counterterms for the top-quark mass (18) and the stop masses (20) are already in a 
general form. The counterterm for the mixing angle (21) is replaced by 



and the counterterm of the A-parameter (23) is rewritten as 
1 



(56) 



5A = — 



U-H. Uh. i^K - ^< ) + (f^*n Ui,, + Ui,, Ui^, ) 6Yi -5m,iA-fi cot (3) . (57) 



Analogous to the top quark/squark sector: 

As in the top quark/squark sector the counterterm for the mixing angle (30) is replaced 



by 



SYr 



2 V 



(5^ 



The dependent counterterms of the 6i-mass (26) and of the A-parameter (31) are rewritten 
as follows: 



SlTlr 



bi: 



■ \-Ul Smf + 2Ur Ur 5Yr + Uf 5mj + U? Sm 

L bi2 b2 bi2 022 hi ti t\2 



5Au 



- "^U-t.^Ui^^SYi + {2nn,6mb - 2mt6mt)] 
1 



(59) 



^^12 Sml + ^5Y, - 5mAA, - // tan/5 - 2^mA 
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(60) 



DR bottom-quark mass 

The A-parameter counterterm (34) is written in the following way 



5A, = — 

nib 



-^ReSr(m, 

U- 022 ^ 02 



2Ur 



bii 



(ReE^;;K,) + ReS^;;(m|)) 



JJ- \ til ill ^ ti' ti2 t22 ^ t; 



m?(Re ES^(m?) + Re Y.'^^im]) + 2 Re Sg^(m?)) ) + ^(2^m, 



2' Ur 



-Ab + ^i tan/3) ( Re E^7(m^) + Re + 2 Re Ksimi) 



div / 



(61) 



avoiding an explicit definition of the mixing angles. The dependent counterterm for the 
mixing angle (35) is replaced by 



U-, 



6Yf^ = j^mbdAb + -^5mh{Ab - /i tan /3 - 2 



U-, 



bi2 



U-, 



Uz 



[6r 
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b22 b22 bii 622 

- U!J< - UlM. + 2Ui^^UiJY-, + 2m,5m,] , 
and the counterterm for the 61-mass (36) by 

5ml =-^\{^-2Ul )(uj 6ml + U? 6ml - 2UiUr6Yi - 2mt6mt) + 6m 

1 U- L t'12'\ ill *1 il2 t2 il2 t22 t ^ ^ ) bx2 ' 

bii 

+ 2U~,^Ui^mb6Ab + (2U-,^^Ui^^{Ab - /x tan/3) + 2(1 - 2UlJmb)6m, 



(62) 



(63) 



DR mixing angle and A5 

The counterterm for the mixing angle (40) is replaced by 

5n = ^(ReSg(m[)+ReS^;;(mij). 



(64) 



The dependent counterterm for the bottom quark mass (41) is rewritten as the following 
combination of counterterms: 



6mh 



6 A, - ^6Y, + ^6ml - ^ (f/|^5m| + f/|/m| - 2U,^^U,JY, 

^bn ^bii ^Sii ^ 



2m+ 5my 



Ah — fi tan P — 2mh 



Uz 



bi2 



Uz 



(65) 



The counterterm for the 61-mass is obtained by inserting the expression (65) for the bottom 
quark mass into the expression (59). 
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On-shell mixing angle and A5 

The renormalised vertex (44) has the following general form 



2M]v sin 

(/i + tan pAb) Unib + l-mb{5Zi^i^ + SZ-^^i^ 



nibtan j36Ab 



(66) 



Using the renormalisation condition (43) the counterterms of the A-parameter and the bot- 
tom quark mass can be derived as 



SAb = i 



Mw sin Ow 



fj. + Ab tan P 



A(0,m?^,m[)+A(0,m|,my) 



2 tan/? 
^ + Ab tan P 
nib tan /? 



Mw sin 9w 



A{0, ml, ml) + A{0, ml, ml 



b2' 



mbif^ + Ab tan/5) ^ ^ ^^^^^ _ ^^^^ -^^iu? . 5m? 



2 tan/? 







1 ■) 

tan/3/ 


+ 2mfe 



til tl 



^4, 

Ur 



(67) 



and 

Smb 



Mw sin ^ 



Letan/?(^,^^^,^^-^,^^f/,^^ 
mfe(yU + Ab tan /?) 



A(0,m[,m?J + A(0, m?^,m|)) 



u-b 

022 



Ur 



2 temp 



is\i.+sZb,-j--Ff^syi ■ —b2 u- 

^feii ^bii ^feii 



+ Ufjml - 2Ui^^UiJYi - 2m,5m, 



ji ( tan P + 



^5m? - 
1 



Uf 5mf 
til ti 



tan P 



2mh 



u 



bi2 



UbJ 



{61 



replacing (45) and (47). The counterterm of the mixing angle (42) is replaced by (58). 
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